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ABSTRACT 

Observations of the energetic neutral atoms (ENAs) of heliospheric origin by IBEX differ from 
expectations based on heliospheric models. It was proposed that the structure of the heliosphere 
may be similar to the ”two-stream” model derived in 1961 by Parker for the case of strong 
interstellar magnetic field. 

Using MHD simulations, we examine possible structure of the heliosphere for a wide range of 
interstellar magnetic field strengths, with different choices of interstellar medium and solar wind 
parameters. For the model heliospheres, we calculate the fluxes of ENAs created in the inner 
heliosheath, and compare with IBEX observations. 

We find that the plasma flow in the model heliospheres for strong interstellar field (~20 /xG) 
has a ”two-stream” structure, which remains visible down to ^5 /xG. The obtained ENA flux 
distribution show the features similar to the ’’split tail” effect observed by IBEX. In our model, 
the main cause of this effect is the two component (fast and slow) solar wind structure. 

Subject headings: Sun: heliosphere — Sun: solar wind — ISM: magnetic fields — magnetohydrodynamics 


1. Introduction 

Energetic neutral atoms (ENAs) created in the 
distant heliosphere by energetic ion neutralization 
provide a means to remotely observe the distant 
regions of the heliosphere. Theoretical models of 
the large scale structure of the heliosphere are im¬ 
portant for understanding and interpreting these 
observations. 

The global models of the stellar wind interac¬ 
tion with the interstellar medium (ISM), leading 
to the formation of the astrospheres (the helio- 
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sphere in the case of the Sun) were firs t introduced 
in the classic work by IParker ( 1961 ). As shown 
by Parker, in the case of a star moving through 
unmagnetized interstellar plasma, the stellar wind 
flow, after passing the termination shock, turns 
ultimately in the direction opposite to the motion 
of the star, forming the ’’tail” (heliotail). This 
structure was indeed obtained in all models of the 
heliosphere based on numerical solutions of the gas 
dynamical or MHD equations. 

Recently, another class of these models was 
included in t he discussion. As again shown by 
Parker! ( 196lll , in the case of a star at rest with re¬ 


spect to the interstellar medium with strong mag¬ 
netic field, the stellar wind may form, instead of 
a single astrotail, two oppositely directed streams 
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or ” drainage plumes”, parallel and antiparallel to 
the magnetic field direction. Although this model 
is not strictly applicable to the case of the Sun 
(which is known to move relative to the local inter¬ 
stellar medium), a structure of this kind was tenta¬ 
tively discussed in the context of the heliospheric 


observations by IBEX and INCA (jMcComas et al 
I 2 OO 9 I : iKrimigis et all 2009l l. 

In a recent study, McComa s et al.l ( 2013) fol- 
lowed on the suggestion of Kivelson fc Jia~(l2013h 


and m ade this idea more explicit. iMcComas et al 


( 201.11 1 discovered that the ENA flux from the 
region of the heliotail in the IBEX data shows 
two regions shifted to opposite sides from the 
downwind direction (anti-apex of the interstellar 
medium flow), and suggested that the splitting 
of the heliotail may be caused by the effect of a 
strong int erstellar magnetic field and subsonic in¬ 
teraction ( McComas et al. 2012ll . The structure 
of the heliosphere would then be somewhere be¬ 
tween the ’’two stream” case (corresponding to 
the extremely strong magnetic field) and the more 
conventional one-tail struc ture, correspon d ing to 
the weak magnetic field. Kivelson fc .lia ( 201, il l 
suggest also that, in the case of strong interstellar 
field, ion heating by reconnection at the heliopause 
may produce a new ENA source. 

Among many numerical models of the helio¬ 
sphere (astrosphere) that take the interstellar 
magnetic field into account, we know of only two 
published models that appr oach close to the ’’two 
stream” extreme field case. Florinski et al.l ( 2004 1 
report a simulation for the case of interstellar 
flow directed along the magnetic field, but they 
found no stationary so lution. For a similar case, 
Pogorelov et al.l ( 201 ill found a solution with the 


heliopause open in both interstellar upwind and 
downwind directions. The case of very strong 
magnetic field oblique to the interstellar flow was, 
however, not considered. As a consequence, there 
are currently no models of the heliosphere that 
co uld be used to examine t he possibility suggested 
by iMcComas et al. ( 2011 1. Such models would 
also be interesting in the case of astrospheres, for 
which different combinations of the magnetic field 
strength and the velocity of the star relative to 
local interstellar medium may be encountered. 

In this work we apply a ID MHD code to obtain 
global models of the heliosphere for a wide range 
of the interstellar field strengths, from 2 /rG to 20 


^G. For the latter value, our models include the 
analogue of the two-stream Parker model when the 
star is at rest. We follow the evolution of the two- 
stream structure when the magnetic field strength 
decreases and its dependence on the neutral hy¬ 
drogen background. 

For these models, we calculate the directional 
distribution of the ENA flux produced inside the 
heliosphere, between the termination shock and 
the heliopause. Our goal is to show how the global 
structure of the heliosphere at different interstellar 
field strengths would be reflected in the ENA ob¬ 
servations. In particular, we find that the plasma 
stream directions correspond to peaks in the ENA 
flux distribution. 


For moderate interstellar field strengths (2 - 
5 fiG) and two-component (slow -|- fast) solar 
wind, our calculations produce the ENA flux de¬ 
pletions reminiscent of the ” split tail” e ffect ob¬ 


served in the IBEX d ata (IMcComas et al 
Schwadron et al.l[2014ll . 


2013; 


We use the following abbreviations: B/s, Vis, 
nis, nn stand for the interstellar medium mag¬ 
netic field strength, plasma speed, proton num¬ 
ber density and neutral hydrogen number den¬ 
sity, respectively. Vsw denotes the solar wind 
speed at the inner boundary and nsw,iAU the so¬ 
lar wind electron density at 1 AU. The ENA flux 
and plasma mass flux directional distributions are 
presented as all-sky maps in solar ecliptic J2000 
coordinates using Mollweide projection with the 
ISM downwind (anti-apex) direction in the center 
and the ISM upwind (apex) direction at both far 
left and far right. 


2. Models of the heliosphere 


Our models of the heliosphere (Figs. [U [2l [I] |4]) 
are defined by time-stationary numerical solutions 
of the (one-fluid) MHD equations for solar wind 
plasma expanding into the magnetized interstellar 
medium. The MHD code used in our calculation 
was shown to be successful in treating the cases of 
different inte rstellar magnetic field strengths and 


orientations (Ratkiewicz et ahlll 


We consider a range of values (2 to 20 fiG) of 
the interstellar magnetic field strength Bjs (Table 
[T]) . Most of the remaining parameters are based on 
observations. The direction of Bis (solar ecliptic 
longitude 221°, latitude 39°) is close to the cen- 
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tre of the IBEX ribbon ( Funsten et alJ[^09l) . The 
direction and magnitude (longitude 79°, latitude 
-4.98°, b/s=23.2 km/s) of the velocity of the inter - 
stellar flow are as given bv lMcComas et ^ ( 2012ll . 
Other interstellar and solar wind parameters (n/s, 
nn, Vswj 'nsw,iAu) are listed in Table [TJ The 
solar wind at the inner boundary is taken either 
spherically symmetric or to consist of two com¬ 
ponents: the slow wind within 36° from the solar 
equator plane, and the fast wind elsewhere. The 
slow and fast c omponents we take to have the same 
ram pressure (jLe Chat. Issautier fc Mever-Vernet 


20121 1. To avoid numerical reconnection, the solar 


magnetic field is neglected (set to zero at the inner 
boundary). 

The neutral gas component is treated as a con¬ 
stant background. The value uh of the order 0.2 
cm“^ i n the inters t ellar m edium is presently fa ¬ 


vored (ll zmodeno^ 


Bzowski et al. ( 2009l l: 


Zank et al.l ( 2013ll l. However, our MHD calcu¬ 
lations assume n//=constant everywhere. If this 
constant is equal to the interstellar value for nn, 
this overestimates the hydrogen density in the in¬ 
ner heliosphere, leading in particular to underes¬ 
timation of the termination shock distance from 
the Sun. For this reason in most of our calcu¬ 
lations we choose n//=0.1 cm“^, in agreement 
with estimations of at the termination shock 
( Bzowski et al. 2009ll rather than 0.2 cm“^. 


Numerical calculations are done on a spheri¬ 
cal (r,0,())) grid with logarithmic spacing for r and 
equal spacing for the angles. The directions of 
the undisturbed interstellar field (Bis) and of the 
inflow velocity of the interstellar matter (Vis) 
together with the position of the Sun define the 
(x,y) plane. The angle 9 is counted from the 
apex direction of the interstellar matter inflow (- 
X axis) and the angle </> from the y axis in the 
(y,z) plane. The numbers of grid points (nr,ne,n^) 
equal (348,90,180) for most calculations. For com¬ 
parison, some calculations were done on smaller 
grids. The calculational domain lies between the 
inner boundary at r=15 AU (in some cases, 30 
AU) and the outer boundary at r=4500 AU. 


3. Models of the energetic ion distribution 
and calculation of the ENA flux 

We are interested in energetic neutral hydrogen 
atoms in the IBEX energy range (^0.7 to ~4.3 


keV) coming from the inner heliosheath (between 
the termination shock and the heliopause). 

The main production mechanism is neutraliza¬ 
tion of the parent ions (energetic protons) by pick¬ 
ing electrons from low energy neutral atoms enter¬ 
ing the heliosphere from the interstellar medium. 
The ENAs from the IBEX ribbon (created pre¬ 
sumably outside the heliosphere) are not consid¬ 
ered. 

Since Voyager 2 observations imply that the 
bulk plasma temperature downstream of the ter¬ 
mination shock is low (~15 eV), we assume that 
most of the parent ions of the ENA derive from 
the pick-up protons created in the solar wind and 
further accelerated near the termination shock. 


Our MHD simulations do not provide the pick¬ 
up ion distribution, so that it must be calculated 
separately. 


First, inside the termination shock we replace 
the constant background model of the neutral hy¬ 
drogen distribution used in the code by the distri- 
bution obtaine d by using a ” hot model” approach 
( ThomasI 1978ll . We take into account the ion¬ 
ization losses, but assume that the gravity is com¬ 
pensated by the radiation pressure. 


The calculated neutral hydrogen distribution is 
used to derive the density of the pick-up protons 
arriving at the termination shock. An example 
of the result is shown in Fig. [S] The density is 
nonuniform, decreasing by a factor of ^2 towards 
the ISM anti-apex direction. 


The distribution function of the energetic pro¬ 
tons at the termination shock we describe by an 
analytical model (see the Appendix). The model 
satisfies the requirement (following from Voyager 2 
observations) that most (0.8) of the solar wind en¬ 
ergy upstream of the shock is transferred to ener¬ 
getic particles instead of heating the bulk plasma. 
Figure m shows the example of the resulting spec¬ 
trum on the opposite sides of the fast/slow wind 
boundary. The pick-up ion density in the fast wind 
region is lower by a factor of ~4-5 than in the slow 
wind. 


The energetic proton flux in the region between 
the shock and the heliopause is then obtained by 
assuming that the particles are convected from the 
shock along plasma streamlines. Losses by neu¬ 
tralization due to charge exchange with ambient 
neutral hydrogen are taken into account. Also in- 
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eluded is the effect of adiabatic acceleration: the 
energy of the energetic particle (in the plasma 
frame) varies along a streamline. 

The ENA flux Jena arriving in the inner solar 
system we calculate by integrating the ENA pro¬ 
duction rate along the lines of sight corresponding 
to the grid directions: 


Jena — / d^Jion^cx^H 


( 1 ) 


where the integral is over the distance s along the 
line of sight (between the termination shock and 
the heliopause), Jena and J^on the fluxes of the 
ENA and of the parent ions, respectively (both 
at the same energy and directed along the line of 
sight), Ucx the charge-exchange cross section at the 
ENA energy (we neglect the speed of the neutral H 
atoms from the background), and tih the number 
density of the neutral H background. 

The characteristic distance for the neutraliza¬ 
tion loss for the energetic protons convected from 
the termination shock is given by V /where V 
is the plasma speed in the inner heliosheath and 
Pcx the rate for charge-exchange between the en¬ 
ergetic proton and the low energy hydrogen atom. 
In the IBEX energy range, /3ca;=0.7-l 10“® s“^. In 
the region of the heliosheath not far from the ter¬ 
mination shock, E~75 km/s for our models with 
ffsw=400 km/s, and E~120 km/s for the models 
with 5 / 5=20 fiG which assume V5w=750 km/s. 
The value of Vjpex is therefore 50-70 AU (400 
km/s solar wind) or 80-110 AU (750 km/h solar 
wind). Most of the production of the ENA from 
the termination shock accelerated protons must 
therefore take place within this distance. Note, 
however, that the fast solar wind effects extend 
this range by increasing V. 


The losses of the ENA on the way to the obser¬ 
vation point are not included in the present calcu- 
l ations, but are s mall for the highe r IBEX energies 
( Bzowski (2008); McComas et al. ( 2012bl l'l 


4. Results: Structure of the model helio¬ 
spheres 

We concentrate on the aspects of the helio¬ 
spheric structure which are most important for 
understanding the production of the ENA in the 
inner heliosheath: ( 1 ) the termination shock and 
the energetic ion distribution in the vicinity of the 


shock, ( 2 ) the structure of the plasma flow down¬ 
stream from the shock (Figures [H [2114|), which de¬ 
termines the transport of the energetic ions from 
the termination shock region into the inner he¬ 
liosheath, and (3) the shape and size of the he¬ 
liopause. 

In Table ffl we collected some results concern¬ 
ing the geometry of the heliosphere for our mod¬ 
els: minimum (rrs.mm) and maximum (rrs.maa:) 
distances to the termination shock, minimum dis¬ 
tance (jHP,min) to the heliopausc, the maximum 
height (’’height”) of the heliosphere counted in 
the z direction (perpendicular to the {Bis,Vis) 
plane), and the width (” width”) of the heliosphere 
along the straight line passing through the point 
of maximum height and perpendicular to the ISM 
inflow direction. 


4.1. Plasma flow structure 


For strong 5 / 5 , the plasma flow between the 
termination shock and the heliopause has a two- 
stream structure. Figures [T] and [ 2 ] show the so¬ 
lar plasma streamlines for 5 / 5=20 /iG starting 
at the termination shock in the symmetry plane 
of the solution (for spherically symmetric solar 
wind). For U/ 5=0 (Figure |T]), the streams are 


respectively parall el and antipa rallel to 5/5 as 
in Parker’s model (jParken 119611) . When Vis 
(Figure 12 ]) the streams appear as bunches of al¬ 
most parallel streamlines and include only a part 
of the flow. The streams are then deflected from 
the ±B /5 directions towards the direction of V /5 
and run approximately parallel to the ’’wings” of 
the heliopause (see next subsection). 

One point to note is the effect of neutral hy¬ 
drogen background on the plasma streams. For 
n//= 0.1 cm“^ (close to the observed value), the 
momentum exchange between the stellar plasma 
and background hydrogen caused by charge- 
exchange interaction deflects and diffuses the 
streams (Fig. |21 second panel). For n//=0.01 
cm“^ and smaller, the two streams are more 
prominent (Fig. |21 first panel) and would be rec¬ 
ognizable even for the weaker field ( 5 / 5=5 ^G). 


A quantitative representation of the two stream 
structure can be obtained by plotting the den¬ 
sity of plasma streamlines. Figure [3] shows the 
directional distribution (projected onto the celes¬ 
tial sphere) of the density of solar plasma stream- 
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lines crossing the Sun-centered sphere of radius 
300 AU for three of our models with spherically 
symmetric solar wind and different values of Bis- 
The streamlines start at the termination shock, 
with their initial points chosen to have the same 
directional distribution as the plasma mass flow. 
The streams appear as two separate density peaks 
near to (but shifted from) the directions parallel 
and antiparallel to the interstellar field. The two- 
stream structure is most prominent at strong Bjs 
(20 /iG) but persists also for Bis=?> /rG and 3 /iG. 


4.2. Shape of the heliopause 

Asymmetric pressure of the interstellar mag¬ 
netic held causes the heliosphere to expand in the 
{Bjs, Vis) plane and contract in the perpendicu¬ 
lar direction (Figure |4]). This effect is well known 
from many numerical simulations. 

For strong Bis (20 fiG and, if ur is small, even 
for 5 ^G) we hnd that the forward part of the 
heliopause in the (i?/s, Vis) plane has the form 
of straight ’’wings” (Fig. [2]). A similar shape 
was predicted using the Newtonian approxima¬ 

tion l^hr. Grzedziels ki fc Rat k iewicz- La iidowskal 
( 198^ . Gzechowski fc Grzedzielskil (|l998 )l. In 
particular, the angle a between the ’’wings” and 
the X axis (the Vis direction) was derived: 


tana = 


tan 7 


1 ± VisIVa cosy’ 


( 2 ) 


where 7 is the angle between the interstellar mag¬ 
netic held and the x axis, and Va is the Alfven 
speed in the interstellar medium. Equation [5] is 
also known from the theory of Alfven wings, which 
appear when a conducting body (like a satellite 
or a planetary magnetos phere) is moving through 
a ma g netized plasma (IDell. Folev fc Ruderman 
(1 1965 1: Neubaueil ( 1980ll : Kivelson fc dial ( 2013 1: 
Saur et 311 1 201311 '). 


In Fig. [21 the directions of the Alfven wings 
are shown by dashed-dotted lines. Although our 
results conhrm the presence of straight ’’wings”, 
the angle between each wing and the x axis is not 
in agreement with Eq. [21 Note that the interstellar 
plasma how in our simulation is subsonic, so that 
the use of Newtonian approximation is not strictly 
justihed. 


4.3. Shape of the termination shock 


For Bis=2Q ^G, and Vis=0 (Figured]), we hnd 
that the shape of the termination shock is very 
close to a Sun-ce ntered sphere , in agreement with 
Parker’s model ( Parker![l961 1. The shock shape 
is weakly elongated along the magnetic held di¬ 
rection, but the difference between the maximum 
and minimum radius of the shock is only 8% of 
the maximum. 

In the Parker model (|Parkei!ll96l[l . the spher¬ 
ical shock appeared as a consequence of assum¬ 
ing incompressible how downstream. In our cal¬ 
culations, the plasma density downstream is not 
constant and the approximate symmetry of the 
shock must have a different explanation. A de¬ 
tailed study of this problem is beyond the scope 
of the present work. 


For Vis we hnd that the ratio between the 
Sun - termination shock distances in the ISM anti¬ 
apex (heliotail) and ISM apex (nose) directions 
goes down when the interstellar held strength in¬ 
creases. The shape of the shock is not much af¬ 
fected by the asymmetry of the solar wind, be¬ 
cause of our assumption that the ram pressures of 
the fast and slow solar wind are equal. 


5. Results: Angular distributions of the 
ENA flux 

The results for directional distribution of the 
ENA hux calculated for the observer in the in¬ 
ner solar system are shown in Eigures |7l [10] and 
[T2l The projections used in our sky maps (Figs. 
[3] [7l [TOl IT^ and IT^ are the same as in Fig. 
6 of iMc Gomas et al.l (12013 ) and in Fig. 10 of 


ISchwadron et al. ( 20140 . that is centered on the 
interstellar downwind direction (white dot in the 
middle). To help with orientation, the position 
of the IBEX ribbon is marked with the line of 
squares, and the interstellar upheld (B/s) and 
downfleld (-B/s) directions are marked by dots. 


5.1. ENA flux distributions for spherically 
symmetric solar wind. 

Figure [7] shows our results for directional distri¬ 
butions of the 4.3 keV ENA flux for Bis=20 //G, 
5 /rG, and 3 ^G assuming spherically symmetric 
solar wind. The distributions are symmetric with 
respect to the (B/s, Vis) plane, which is shown 
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as a thick black line. For the case of i?/s=20 
and Bis=5 fxG, the ENA flux has two peaks at 
directions approximately parallel and antiparallel 
to the interstellar field, corresponding to the two 
streams of the solar plasma. For Bjs=S /J.G, these 
peaks almost completely disappear. 

The reason why (for Bjs=20 and 5^G) the 
two solar plasma streams are associated with high 
ENA flux intensity can be explained as follows (see 
Figs. [51E]). The lines-of-sight directed along the 
streams are close to parallel to the local plasma 
flow (Fig. E]), so that Eq. IBBI is applicable. Since 
PcxL/V ^ 1 in the stream regions, it follows that 
Jena ~ {V/v)Jo. This relation agrees well with 
our results for the case of Bis=5 /iG without adi¬ 
abatic acceleration. 

The same argument can be applied to the lines- 
of-sight in the heliotail. We find that the ENA flux 
from these directions is lower than for the streams 
because V and Jq are lower for the case of the 
heliotail. 

For Bis=5 and 3 /rG, a broad maximnm of the 
ENA flux appears in the ’’nose” (ISM apex) re¬ 
gion (Fig. [3 panel 2 and 3). This is caused by 
the adiabatic acceleration of the energetic protons 
downstream from the shock and disappears if the 
adiabatic acceleration is switched off (Fig. (3 panel 
4 and 5). 

Eq. IB4I provides the estimation for the ENA 
flux from the ’’nose” direction in the absence of 
adiabatic acceleration. For Bis=5 fiG, the pa¬ 
rameter /3ca;T/Vb~l/2, implying that Jena ~ 
(l/3)(Vb/'c) Jq. This agrees well with our numer¬ 
ical result when the adiabatic acceleration is dis¬ 
regarded (Jo=110 cm“^s“^sr“^keV“^, Jena=^-8 
cm“^s“^sr“^keV“^). 

5.2. ENA flux distributions for asymmet¬ 
ric solar wind: The ’’split tail” for 
high energy ENA flux 

The first three years of IBEX data refer to 
the extended solar minimum period, when the so¬ 
lar wind asymmetry (two-component slow -|- fast 
structure) was prominent. The results of our cal¬ 
culations for the case of asymmetric solar wind are 
collected in Figures [13 [H] and (TOl 

For very strong field (i?/s=15 /rG, Fig. [TT]), the 
peaks from the two plasma streams are prominent. 

For Bis=5 ^J‘G and less, the directional distri¬ 


bution of the ENA flux following from our models 
reflects primarily the two-component solar wind 
structure. Upstream from the shock, the slow 
wind is in our models restricted to the region be¬ 
tween ±36° solar latitude (red lines in Figs. [TUI and 
[T2l) . Figure [To] shows that, at low ENA energy (1.1 
keV) this region corresponds to high ENA inten¬ 
sity belt. At high ENA energy (4.3 keV. Fig. [T2ll 
there is a similar belt, but of low ENA intensity. 


This structure can be understood by observ¬ 
ing that the ENA coming from directions within 
the ±36° solar latitude derive mostly from the en¬ 
ergetic protons carried by the shocked slow solar 
wind plasma. The density of ^^1 keV energetic 
protons is much higher in the shocked slow wind 
than in the shocked fast wind (Fig. [6]). This ex¬ 
plains the high ENA intensity belt at low (~1 keV) 
energy. On the other hand, the energetic proton 
flux at ~4.3 keV is approximately the same for the 
slow and fast solar wind plasma, because the high 
average energy of the energetic protons in the fast 
wind is compensated by the high density of the 
slow wind (Fig. [5]). The reason why the 4.3 keV 
ENA flux is higher in the region outside the slow 
wind belt (Fig. [TU|) is that the convection time 
from the termination shock (and consequently the 
neutralization loss) for the energetic protons in the 
fast wind plasma is lower than in the slow wind. 


For Bis 2-5 fj,G (Fig. [T2|), we find that the 
low intensity belt for 4.3 keV ENA narrows down 
near the ISM downwind direction. This structure 
is very much alike to th e ’’spli t tail” observed by 


IBEX ([McGomas et al.l (|2013ri . ISchwadron et al 


(120141 11. The low ENA flux region around the 
ISM downwind direction is almost split into two 
parts (the ’’po r t” an d the ’’starboard” lobes: 
McGomas et al.l (201^). We find that this ef¬ 


fect is caused by the streams of fast wind plasma 
entering the heliotail region. Figure [T2| (second 
panel) shows the outline of the heliopause (thick 
black oval) at the distance 250 AU from the Sun, 
together with the outlines of the two regions in 
the heliotail filled by the plasma originating from 
the fast solar wind (thin black ovals). The flux of 
~4 keV ENA from these fast wind plasma streams 
is higher than from the slow wind plasma region 
between them. Since the fast wind streams at 
large distances move to lower heliolatitudes, the 
two high ENA flux regions associated with them 
intrude into the ±36° heliolatitude belt and partly 


6 













split the low ENA flux region. 

The low 4.3 keV ENA flux region obtained from 
our calculations (Fig. [12]) is tilted relative to eclip¬ 
tic. This tilt is not dependent on Bjs value (see 
Fig. (HI). On the other hand, it is close to the 
7.25° tilt of the solar equator (and of the slow/fast 
solar wind boundary: see the red lines in the fig- 
ures). A similar tilt a ppears in the IBEX data 
( McComas et al.l ( 2013 ). Fig. 6 ). It may be con¬ 
cluded that the main cause of the ’’split tail” struc¬ 
ture is the distribution of the slow and fast solar 
wind streams. 

Our results concerning the ” split tail” structure 
remain valid when the effect of adiabatic acceler¬ 
ation is switched off (Fig. [T^l two lower panels). 
This is important, because we expect that adia¬ 
batic acceleration is overestimated by our models. 
On the other hand, adiabatic acceleration of the 
energetic protons affects the production of ENA 
from the region near the apex of the ISM (the 
nose region: Fig. dl two upper panels). 

5.3. Comparison with IBEX 

In addition to the ” split tail” structure, our re¬ 
sults show other points of qualitative similarity 
with IBEX. At low energy, the ENA flux in our 
calculations have a peak near the tail direction 
(Fig. [T 0 |l which, at higher energy, evolves into a 
depletion (Fig. IT^ . A similar behaviour can be 


(I2014D . 


seen in the IBEX data (iSchwadron et al. 

Fig. 10). 

Compared to the ENA fluxes observed by 
IBEX, our calculations produce lower ENA in¬ 
tensity. The magnitude of the ENA flux at 4.3 
keV obtained in our models is lower than the non¬ 


of 2 (Fig. 6 in 

McComas et al.l ( 

20131). Fig. in 

in Schwadron et al.l (2014!)). For 

^1.1 keV this 


factor is about 4. 

In Figure [T4| we show our results for the slope of 
the ENA flux spectrum for the case of 5 fiG held, 
asymmetric solar wind. We plot the values of 7 , 
dehnedasy = - \og{ Jena{Ei)/Jena{E 2))/\og{Ei/E2) 
for {Ei,E 2 )=il.l keV,1.7 keV) and (2.7 keV, 

4.3 keV). The values of the spectral index ob¬ 
tained in our model ar e lower t han observed by 


IBEX ( McComas et al. ( 20131) . Schwadron et al 


(|2014D ). suggesting that our model underestimates 
the ENA hux at low energy compared to the high 


energy hux. On the other hand, we hnd that the 
high 7 region is similar to the ” split tail” structure 
in the ENA hux di stribution, in agreeme nt with 
IBEX observations (iMcComas et al.l ( 20131) . Fig.5; 


Schwadron et al.l (|2014l) . Fig. 11, top panel). 


6. Summary and conclusions 

We present a set of 3D time-stationary models 
of the heliosphere based on numerical MHD solu¬ 
tions for a wide range of interstellar magnetic held 
strength ( 2-20 /rC). 

For the ISMF of 20 fiG and the star at rest, 
our result is analogous to a well known Parker so¬ 
lution, with a spherical termination shock, and the 
astrosphere elongated along the Bjs direction. 

For strong interstellar held and Vis 0, we 
show that the plasma how inside the heliosphere is 
concentrated in two streams directed close (though 
not exactly parallel or antiparallel) to the inter¬ 
stellar held direction. As a result, the forward 
part of the heliopause forms straight ’’wings”. A 
similar shape was predicted using the Newtonian 
approximation. 

When the held strength decreases to realis¬ 
tic values (2-5 /rC), the two-stream structure be¬ 
comes less prominent. This structure depends on 
the interstellar neutral hydrogen background: the 
weaker nn, the more prominent the streams. 

For the simulated heliospheres we calculated 
the distribution of the energetic protons (originat¬ 
ing from the pick-up protons from the solar wind) 
and the energetic neutral atom huxes produced 
by neutralization of these protons. We assumed a 
simple model of the energetic proton distribution 
at the termination shock. The results were com¬ 
pared with the ENA observations by IBEX. Our 
simulations are restricted to the ENA created in¬ 
side the heliosphere, so that the ’’ribbon” contri¬ 
bution must be excluded from the IBEX data. 

We find that the structure of the ENA flux 
distribution at higher IBEX energy (1.7 keV and 
higher) in the downwind hemisphere is similar to 
the ’’split tail” structure observed by IBEX (two 
low ENA flux ” lobes” shifted relative to the down¬ 
wind direction). This result persists for different 
values (2-5 /rG) of the interstellar magnetic field 
strength. The explanation suggested by our model 
is that this effect follows from the two component 
(fast -|- slow) solar wind structure, prevalent near 
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the low activity parts of the solar cycle. The tilt 
of the observed ENA flux structure is explained 
by the tilt of the slow/fast solar wind boundary 
relative to the ecliptic plane. 


For spherically symmetric solar wind, as ob¬ 
served near solar maxima, the ENA distribution 
in our time-stationary model has a different form. 
In reality, because of long time delay (average ^3 
years) before neutralization of the parent proton, 
we expect that the ENA distribution near solar 
maximum would still show some effects of the two- 
component solar wind. 

For very strong Bjs, the ENA distributions 
have two prominent peaks corresponding to two 
plasma streams. 

There are also peaks of the ENA flux that ap¬ 
pear in the regions near the nose of the helio¬ 
sphere. However, these effects disappear when 
adiabatic acceleration of energetic protons is ne¬ 
glected. Adiabatic acceleration depends on details 
of plasma flow and density distribution and is sen¬ 
sitive to numerical effects. The ’’split tail” struc¬ 
ture and the two stream-related ENA peaks are 
not dependent on adiabatic acceleration. 


Finally, we note that the recen t time-dependent 
simulation ( Zirnstein et al. 2015[l found the ENA 
flux distribution from the heliosphere with little 
periodic change over the solar cycle. Their ap¬ 
proach was different from ours. In particular, 
the pick-up protons were not treated separately, 
but assumed to form a fixed fraction of the bulk 
plasma. 


We thank M. Strumik for many helpful dis¬ 
cussions. Work in the U.S. was supported by 
the IBEX mission as part of NASA’s Explorer 
Program. A.C. was supported by the grant 
2012/06/M/ST9/00455 from the Polish National 
Science Center. 







A. Energetic ion distribution 


The method is similar to that used in Czechowski et al. (jCzechowski. Hilchenbach fc Hsieh (j2012li i. The 
velocity distribution of energetic protons f{f,v) is assumed to be isotropic in plasma frame. It is obtained 
by solving the time-stationary transport equation 


[-v-^-Pc. + Uv-v)j^]f = o 

3 o log u 


(Al) 


where V is the plasma velocity, v the particle velocity and f3cx is the rate for charge exchange between the 
energetic protons and low energy neutral hydrogen. The equation includes the effects of convection, adiabatic 
energy changes and neutralization loss. We assume that the proton energy is low enough to neglect spatial 
diffusion. 

Equation lAll together with the boundary condition at the termination shock determine the f{r,v) along 
a plasma streamline r = r(s) parametrized by the length s: 


f{f{s),v{s)) = fshock{ro,vo)exp 


/o ^(r{sO)J 


(A2) 


where r{s) and v{s) satisfy the equations VAr/As = V{s) and Edlogw/ds = —(1/3)V • V, with the initial 
conditions at the termination shock r(0)=ro, z;(0)=?;oj respectively. If the plasma mass flow is conserved, 
the equation for v(s) implies u(s)p(s)“^/^=const. Since determination of V • E from numerical output is 
imprecise, we use the above relation to determine v(s). 

We use a simple analytical model for the energetic proton distribution at the termination shock /shock (r, v) 
in the plasma frame 

/shock(r, 'c) = npui(r)-FK(v, w) -I- nswif^Fciv, vt) (A3) 


The first term is the distribution of the accelerated pick-up protons (number density npui), and the second 
describes the protons from the shock-heated bulk solar wind (number density nsw)- Fk{v,w) is the kappa 
function 


F^{v,w) 


r(At-H) 1 

r(K — 1 / 2 ) 


2 \ 


1 + 




(A4) 


We choose k= 1.65. The ’’thermal speed” parameter w we determine from the requirement that the total 
pick-up proton energy should be equal to 0.8 of the energy of the solar wind upstream from the shock, as 
observed by Voyager 2 ( Richardson et al.1l2008 ). The bulk solar wind protons are described by the Maxwellian 
distribution Fc{v, vt) = exp —{v/vt)‘^/t^^^‘^v^ with the thermal speed vt equal to one half of the solar wind 
speed downstream. The number density nsw is taken from the MHD model. 


B. ENA flux estimations 

Assume that the plasma flow is incompressible, so that Eq. lAll becomes 


V • V/ = -/3c./ 


(Bl) 


The same equation holds for the ion flux Jion- 

Consider first the ENA flux from the direction towards the nose of the heliosphere, along the stagnation 
line. Let the plasma velocity decrease linearly towards the heliopause: V = Vb(l — z/L) where z is the 
distance from the shock along the stagnation line and L the distance from the shock to the heliopause. 


The equation for the ion flux is 


dJir 


Pc. 


dz 


V^{l-zlL) 


die 


(B2) 
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with the solution 


J^on[z) = Jo (1 - 


The ENA flux is then given by (see Eq. [T]) 


Jena 


PcxL/Vq Eq 
1 + V 


(B3) 


(B4) 


where v is the particle speed and we have used /3cx = o’cxvnn- 

Consider next the line-of-sight within one of the plasma streams assuming that the plasma speed V is 
constant and parallel to the line-of-sight. The solution for the ion flux is then 


Jion{z) = Jo exp(-2;/3ca;/14) 


(B5) 


and for the ENA flux 

•/ejva =-^Jo -exp (B6) 

When the distance to the heliopause along the line-of-sight is large compared to Vjficx: the ENA flux becomes 
Jena = iV/v)Jo. 


10 




REFERENCES 

Burlaga, L.F., Ness, N.F., 2014, ApJ, 784, 146 

Bzowski, M., 2008, Astronomy & Astrophysics, 
488,1057 

Bzowski, M., Mobius, E., Tarnopolski, S., Izmode- 
nov, V., Gloeckler, G., 2009, Space Science Re- 
vievs, 143, 177 

Bzowski, M., Kubiak, M.A., Mobius, E., Bochsler, 
P., Leonard, T., Heirtzler, D., Kucharek, 
H., Sokol, J.M., Hlond, M., Crew, G.B., 
Schwadron, N.A., Fuselier, S.A., McComas, 
D.J., 2012, ApJS, 198, 12 

Czechowski, A., Grzedzielski, S., 1998, Geo- 
phys. Res. Lett., 25, 1855 

Czechowski, A., Hilchenbach, M., Hsieh, K.C., 
2012 , Astronomy & Astrophysics, 541, 14 

Drell, S.D., Foley, H.M., Ruderman, M.A., 1965, 
J. Geophys. Res., 70, 3131 

Fahr, H.J., Grzedzielski, S., Ratkiewicz- 
Landowska, R., 1988, Ann. Geophys., 6, 
337 

Florinski, V., Pogorelov, N.V., Zank, G.P., Wood, 
B.E., Cox, D.P., 2004, ApJ, 604, 700 

Funsten, H. O., Allegrini, F., Crew, G. B., DeMa- 
jistre, R., Frisch, P. C., et ah, 2009, Science, 
326, 964 

Gruntman, M., Roelof, E.C., Mite hell, D.G., 
Fahr, H.J., Funsten, H.O., McComas, D.J., 
2001, J. Geophys. Res., 106, 15767 

Grygorezuk, J., Czechowski, A., Grzedzielski, S., 
2014, ApJ, 789, L43 

Heerikhuisen, J., Pogorelov, N., Zank, G.P., Crew, 
G.B., Frisch, P.C., Funsten, H.O, Janzen, P.H., 
McComas, D.J., Reisenfeld, D.B., Schwadron, 
N.A., 2010, ApJ, 708, L126 

Heerikhuisen, J., Zirnstein, E.J., Funsten, H.O, 
Pogorelov, N.V., Zank, G.P., 2014, ApJ, 784:73, 
doi:10.1088/004-637X/784/l/73 

Izmodenov, V.V., 2009, Space Sci. Revs., 143, 139 

Kivelson, M.G., Jia, X., 2013, J. Geophys. Res., 0, 
0 (submitted) 


Krimigis, S.M., Mitchell, D.G., Roelof, E.C., 
Hsieh, K.C., McComas, D.J., 2009, Science, 
326, 971 

Le Chat, G., Issautier, K., Meyer-Vernet, N., 2012, 
Solar Phys., 279, 197 

McComas, D.J., F. Allegrini, P. Bochsler, M. 
Bzowski, E.R. Christian, G.B. Crew, R. De- 
Majistre, H. Eahr, H. Eichtner, P.C. Frisch, 
H.O. Funsten, S.A. Fuselier, G. Gloeckler, M. 
Gruntman, J. Heerikhuisen, V. Izmodenov, P. 
Janzen, P. Knappenberger, S. Krimigis, H. 
Kucharek, M. Lee, G. Livadiotis, S. Livi, R.J. 
MacDowall, D. Mitchell, E. Mobius, T. Moore, 
N.V. Pogorelov, D. Reisenfeld, E. Roelof, L. 
Saul, N.A. Schwadron, P.W. Valek, R. Vander- 
spek, P. Wurz, G.P. Zank, 2009, Science, 326, 
959 

McComas, D.J., Alexashov, D., Bzowski, M., 
Fahr, H., Heerikhuisen, J., Izmodenov, V., Lee, 

M. A., Mobius, E., Pogorelov, N., Schwadron, 

N. A., Zank, G.P., 2012, Science, 336, 1291 

McComas, D.J., Dayeh, M.A., Allegrini, E., 
Bzowski, M., DeMajistre, R., Fujiki, K., Fun¬ 
sten, H.O., Fuselier, S.A., Gruntman, M., 
Janzen, P.H., Kubiak, M.A., Kucharek, H., 
Livadiotis, G., Mobius, E., Reisenfeld, D.B., 
Reno, M., Schwadron, N.A., Sokol, J.M., Toku- 
maru, M., 2012, ApJS, 203, 1 

McComas, D.J., Dayeh, M.A., Eunsten, H.O., Li¬ 
vadiotis, G., Schwadron, N.A., 2013,ApJ, 771, 
77 

Mobius, E., Bochsler, P., Bzowski, M., Heirtzler, 
D., Kubiak, M.A., Kucharek, H., Leonard, T., 
Lee, M.A., Leonard, T., Schwadron, N.A., Wu, 
X., Fuselier, S.A., Crew, G., McComas, D.J., 
Petersen, L., Saul, L., Valovcin, D., Vander- 
spek, R., Wurz, R, 2012, ApJS, 198, 11 

Neubauer, F.M., 1980, J. Geophys. Res., 85, 1171 

Parker, E.N., 1961, ApJ, 134, 20 

Pogorelov, N.V., heerikhuisen, J., Zank, G.P., 
Borovikov, S.N., Frisch, P.C., McComas, D.J., 
2011, ApJ, 742, 104, doi: 10.1088/0004- 
637X/742/2/104 

Richardson, J.D., Kasper, J.C., Belcher, J.W., 
Lazarus, A.J., 2008, Nature, 454, 63-66 


11 



Ratkiewicz, R., Barnes, A., Molvik, G. 

A.,Spreiter, J. R., Stahara, S. S., Vinokur, M., 
and Venkatesvaran, S. A., 1998, A&A, 335, 
363-369 

Saur, J., Grambusch, T., Duling, S., Neubauer, 
F.M., Simon, S., 2013, Astr. Astrophys., 552, 
A119 

Schwadron, N.A., Bzowski, M., Crew, G.B., 
Gruntman, M., Fahr, H., Fichtner, H., 
P.C. Frisch, H.O. Funsten, S. Fuselier, J. 
Heerikhuisen, V. Izmodenov, H. Kucharek, M. 
Lee, G. Livadiotis, McComas, D.J., E. Mbbius, 
T. Moore, J. Mukherjee, N.V. Pogorelov, C. 
Prested, D. Reisenfeld, E. Roelof, G.P. Zank, 
2009, Science, 326, 966 

Schwadron, N.A., Allegrini, E., Bzowski, M., 
Christian, E.R., Crew, G.B., Dayeh, M., DeMa- 
jistre, R., Erisch, P., Eunsten, H.O., Euselier, 

5. A., Goodrich, K., Gruntman, M., Janzen, P., 
Kucharek, H., Livadiotis, G., McComas, D.J., 
Mobius, E., Prested, C., Reisenfeld, D., Reno, 
M., Roelof, E., Siegel, J., Vanderspek, R., 2011, 
ApJ, 731, 56 

Schwadron, N.A., Mobius, E., Euselier, S.A., 
McComas, D.J., Funsten, H.O., Janzen, P., 
Reisenfeld, D., Kucharek, H., Lee, M.A., 
Fairchild, K., Allegrini, F., Dayeh, M., Liva¬ 
diotis, G., Reno, M., Bzowski, M., Sokol, J.M., 
Kubiak, M.A., Christian, E.R., DeMajistre, R., 
Frisch, P., Galli, A., Wurz, P., Gruntman, M., 
2014, ApJS, 215, 13 

Stone, E.C., Cummings, A.C., McDonald, F.B., 
Heikkila, B.C., Lai, N., Webber, W.R., 2008, 
Nature, 454, 71-74 

Thomas, G.E., 1978, Ann. Rev. Earth Planet. Sci. 

6 , 173-204 

Wood, B.E., Izmodenov, V.L, Alexashov, D.B., 
Redfield, S., Edelman, E., 2014, ApJ, 780, 108 

Zank, G.P., Heerikhuisen, J., Wood, B.E., 
Pogorelov, N.V., Zirnstein, E., McComas, D.J., 
2013, ApJ, 763, 20 

Zirnstein, E.J., Heerikhuisen, J., Pogorelov, N.V., 
McComas, D.J., Dayeh, M.A., 2015, ApJ, 804, 
5 


This 2-column preprint was prepared with the AAS lAT^ 


macros v5.2. 


12 




B = 20/tG V=0 n = 0.04 nH = 0.1 

600 

' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 

400 

^ _ 

200 

■ ■ 

S ° 

- 


: : 

-200 


-400 

: ^ : 


. : 

-600 

- 


_1 ... 1 ... 1 ... 1 ... 1 ... 1 . • 


-600 -400 -200 0 200 400 600 

x(AU) 


Fig. 1.— Flow lines of the solar plasma for the 
MHD solution corresponding to the hypothetical 
heliosphere with the Sun at rest (yfs=0) relative 
to the interstellar medium with plasma density 
n/s=0.04 cm“^, neutral hydrogen density nij=0.1 
cm“^ and very strong interstellar magnetic field 
(5/5=20 ^G). The solar wind is spherically sym¬ 
metric with Vsvv=750 km/s, n^wA 4r/=4.2 cm~^. 
The solution is similar to the Parker model (|Parker 


19611) . The heliopause and the termination shock 


are shown by thick lines. 


B = 20/xG V=23.2 n = 0.04 n„=0.01 




X (AU) 

Fig. 2.— Plasma flow lines in the {x,y)={Bis,Vis) 
plane obtained from MHD simulations for 5/s=20 
/iG, V/s=23.2 km/s and the solar wind with 
ksw=750 km/s, nsrv,iA(7=4.2 cm“^, The first fig¬ 
ure corresponds to n//=0.01 cm“^ and n/s=0.04 
cm“^, and the second to n//=0.1 cm“^ and 
n/s=0.06 cm“^. The solar wind is spherically 
symmetric, so that the (x,y) plane is the symme¬ 
try plane. The heliopause is shown by the thick 
line. The dotted line shows the 5/s direction. The 
dashed lines show the Alfven wings (see Eq. [2]). 
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Fig. 3.— Directional distribution of the solar 
plasma mass flux (obtained from streamline den¬ 
sity) in units 10^ rup cm“^ s“^ at a distance 300 
AU from the Sun for the cases of Bjs=20 (the 
first panel) 5 fiG and 3 fiG. The values of D/s=23.2 
km/s, n//=0.1 cm“^, n/s=0.06 cm“^ are the same 
for all cases. The solar wind is spherically sym¬ 
metric, with Vsw=750 km/s, nsw,i7i;7=4.2 cm“^ 
(upper panel) and Vsw=400 km/s, nsw,1/1(7=5.55 
cm“^ for the remaining cases. The projection on 
the celestial sphere is the same as used for the 
ENA distributions (Figs. IT^ ITOl and iTl) . It is cen¬ 
tered on the anti-apex direction of the ISM flow 
(white circle). The interstellar field and anti-field 
directions are marked by black circles. 


Fig. 4.— Plasma velocity distributions in the (x,y) 
plane (first panel) and the (x,z) plane (second 
panel) for Bis=h ^G, n//=0.1 cm“^, D/s=23.2 
km/s, n/s=0.06 cm“^. The solar wind is 
asymmetric (fast/slow) with Dsv^=750/400 km/s, 
7isw,i/ic/=l-58/5.55 cm“^. The outlines of the 
heliopause and of the termination shock are also 
shown. Note the flattening of the heliosphere by 
the asymmetric pressure of the interstellar mag¬ 
netic field. 


14 


400/750 
































Fig. 5.— The pick-up proton density at the termi¬ 
nation shock as a function of angle 9 (counted from 
the inflow direction of the interstellar medium) 
for two choices of the angle (j) corresponding to 
maximum (solid line) and minimum (dashed line) 
densities for intermediate 9. The cases shown 
are: Bis=5 /xG, symmetric solar wind (Gsu^=400 
km/s, nsvix,iAf7=5.55 cm“^, the upper panel) and 
Bis=5 /xG, asymmetric (fast/slow) solar wind 
(1(5^=750/400 km/s, nsvy,iAC/=1.58/5.55 cm-^, 
lower panel). The low pick-up proton density re¬ 
gion in the lower figure (dashed line) corresponds 
to the fast wind. The distributions have a similar 
form for other values of Bis, provided that Vis is 
the same (23.2 km/s). 



0.01 0.10 1.00 10.00 
ENERGY (keV) 


Fig. 6.— The proton flux as a function of en¬ 
ergy at two selected locations on the termination 
shock, corresponding to slow (solid line) and fast 
(dashed line) solar wind, respectively. The case 
of Bis=5 /xG, asymmetric (fast/slow) solar wind 
(1/5^=750/400 km/s, nsw.iAC/=1.58/5.55 cm-^ 
is illustrated. The flux of high energy protons up 
to ~10 keV is higher in the slow solar wind region 
because of high slow solar wind density. 
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Fig. 8.— Streamlines ending at selected lines-of- 
sight within 250 AU from the Sun. Two of the 
lines-of-sight correspond to peaks of the ENA flux 
associated with two plasma streams. In addition, 
one line-of-sight is in the forward part of the he¬ 
liosphere (the ’’nose” region) and one is close to 
the ISM anti-apex direction (heliotail). The case 
illustrated corresponds to the 1st panel in Fig. [7] 
{Bjs=20 ^G, symmetric solar wind). 


Fig. 7.— ENA flux distribution (4.3 keV) in units 
(cm^ s sr keV)“^ for spherically symmetric solar 
wind. The cases shown (from left to right and from 
top to bottom) are: i3/s=20 fiG, 5 fiG, 3 fiG, 5 
fiG and 3 fiG, the two last cases corresponding 
to adiabatic acceleration switched off. The ENA 
flux for Bis=20 fiG and 5 /xG has two peaks corre¬ 
sponding to two streams of the plasma flow. These 
streams effectively disappear for 3 /rG. The ENA 
flux from the ISM apex region is affected by adia- 16 
batic acceleration. The values of V/s=23.2 km/s, 
n 7 s= 0.06 cm“^, 7177 = 0 .! cm“^, U5vy=400 km/s, 
nsw,iAU=5-55 cm“^ are the same for all figures. 















5 /xG sym 



Fig. 9.— Profiles of the ENA production rate 
(4.3 keV) along 5 different directions correspond¬ 
ing to: ENA flux maximum in the south hemi¬ 
sphere; ENA flux maximum in the north hemi¬ 
sphere; vicinity of the ISM anti-apex; the middle- 
point between the two maxima of the flux. The 
case illustrated corresponds to the 4th panel in 
Fig. [7] (5/5=5 fiG, symmetric solar wind, adia¬ 
batic acceleration switched off). 




Fig. 10.— ENA flux distribution (1.1 keV) in 
units (cm^ s sr keV)“^ for Bjs=5 //G (first panel) 
and Bis=3 fiG (second panel), for the case of 
asymmetric (slow-|-fast) solar wind, with slow 
wind contained within ±36° from the equator. 
The values of V/s=23.2 km/s, n/s=0.06 cm“^, 
n//=0.1 cm“^, Fsrv=750/400 km/s (fast/slow), 
»^Svy,i/u/=l-58/5.55 cm“^ (fast/slow) are the 
same for both cases. 



Fig. 11.— ENA flux distribution (4.3 keV) in units 
(cm^ s sr keV)“^ for 5/s=15 //G. 
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Fig. 12.— ENA flux distribution (4.3 keV) in 
units (cm^ s sr keV)“^ for Bjs 5 and 3 fiG. 
The solar wind is asymmetric (fast /slow) with the 
slow wind contained within ±36° from the equa¬ 
tor. The thin red lines are the projections of 
the solar equator and of the boundaries between 
the slow and fast solar wind. The second fig¬ 
ure includes the outline of the heliopause and of 
the regions filled by the shocked fast solar wind 
plasma at the distance 250 AU from the Sun. The 
last two figures show the results {Bis= 5 and 
3 /xG) with adiabatic acceleration switched off. 
The values of 14/5=23.2 km/s, n/5=0.06 cm“^, 
n//=0.1 cm“^, t45M/=750/400 km/s (fast/slow), 
fT-Sw.iAC/=l-58/5.55cm“^ (fast/slow) are the same 
for all figures. 
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Fig. 13.— Profiles of the ENA production rate 
(4.3 keV) along 4 different directions correspond¬ 
ing to: the port lobe; the starboard lobe; the 
’’pinch” from the north, and the ’’pinch” from the 
south. The case illustrated corresponds to the low¬ 
est panel in Fig. [T2](i?/s=3 //G, asymmetric solar 
wind, adiabatic acceleration switched off). 



'.0 


Fig. 14.— Spectral index distributions in the 
1.1-1.7 keV (first panel) and in the 2.7-4.3 keV 
(second panel) energy ranges for numerical so¬ 
lution with 73/5=5 fiG, asymmetric (fast±slow) 
solar wind, n//=0.1 cm“^, 14/5=23.2 km/s, 

71/5=0.06 cm“^, Vsu/=750/400 km/s (fast/slow), 
^SW,i/u/=l-58/5.55cm“^ (fast/slow). The pro- 
j ection is the same as in Figs. 5 and 6 of 
( McGomas et al.ll2013ll l. 
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Table 1: Parameters and results of MHD calcula¬ 
tions 


Bis 

{^^G) 

nis 

(cm“^) 

Uh 

cm“^ 

Vis 

(km/s) 

Vsw 
(km/s) 
fast /slow 

nsw,iAu 
(cm“^) 
fast / slow 

rTS.min 

(AU) 

^TS,max 

(AU) 

rHP.min 

(AU) 

height 

(AU) 

width 

(AU) 

20 

0.04 

0.1 

0 

750 

4.2 

36 

40 

95 



20 

0.06 

0.1 

23.2 

750 

4.2 

36 

43 

57 

140 

749 

20 

0.06 

0.01 

23.2 

750 

4.2 

35 

41 

65 

154 

1840 

15 

0.06 

0.1 

23.2 

400 

5.55 

31 

38 

47 

113 

444 

15 

0.06 

0.1 

23.2 

750/400 

1.58/5.55 

31 

39 

45 

130 

483 

5 

0.06 

0.1 

23.2 

400 

5.55 

63 

91 

97 

282 

622 

5 

0.06 

0.001 

23.2 

400 

5.55 

75 

109 

109 

317 

1827 

5 

0.06 

0.1 

23.2 

750/400 

1.58/5.55 

62 

92 

93 

315 

693 

3 

0.06 


23.2 

400 

5.55 

75 

122 

115 

384 

643 

3 

0.06 

0.1 

23.2 

750/400 

1.58/5.55 

74 

120 

111 

424 

719 

2 

0.06 

0.1 

23.2 

400 

5.55 

81 

140 

119 

469 

654 

2 

0.06 

0.1 

23.2 

750/400 

1.58/5.55 

76 

135 

111 

536 

720 
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